Sea ice assimilation into a coupled ocean–sea ice model using its adjoint

Satellite sea ice concentrations (SICs), together with several ocean parameters, are assimilated into a regional Arctic coupled ocean–sea ice model covering the period of 2000–2008 using the adjoint method. There is substantial improvement in the representation of the SIC spatial distribution, in particular with respect to the position of the ice edge and to the concentrations in the central parts of the Arctic Ocean during summer months. Seasonal cycles of total Arctic sea ice area show an overall improvement. During summer months, values of sea ice extent (SIE) integrated over the model domain become underestimated compared to observations, but absolute differences of mean SIE to the data are reduced in nearly all months and years. Along with the SICs, the sea ice thickness fields also become closer to observations, providing added value by the assimilation. Very sparse ocean data in the Arctic, corresponding to a very small contribution to the cost function, prevent sizable improvements of assimilated ocean variables, with the exception of the sea surface temperature

Simple algorithms to compute meridional overturning and barotropic streamfunctions on unstructured meshes

Abstract. Computation of barotropic and meridional overturning streamfunctions for models formulated on unstructured meshes is commonly preceded by interpolation to a regular mesh. This operation destroys the original conservation, which can be then artificially imposed to make the computation possible. An elementary method is proposed that avoids interpolation and preserves conservation in a strict model sense. The method is described as applied to the discretization of the Finite volumE Sea ice – Ocean Model (FESOM2) on triangular meshes. It, however, is generalizable to colocated vertex-based discretization on triangular meshes and to both triangular and hexagonal C-grid discretizations. </jats:p

On the Damping Time Scale of EVP Sea Ice Dynamics

We propose to make the damping time scale, which governs the decay of pseudo-elastic waves in the Elastic Viscous Plastic (EVP) sea-ice solvers, independent of the external time step and large enough to warrant numerical stability for a moderate number of internal time steps. A necessary condition is that the forcing on sea ice varies slowly on the damping time scale, in which case an EVP solution may still approach a Viscous Plastic one, but on a time scale longer than a single external time step. In this case, the EVP method becomes very close to the recently proposed modified EVP (mEVP) method in terms of stability and simulated behavior. In a simple test case dealing with sea ice breaking under the forcing of a moving cyclone, the EVP method with an enlarged damping time scale can simulate linear kinematic features which are very similar to those from the traditional EVP implementation, although a much smaller number of internal time steps is used. There is more difference in sea-ice thickness and linear kinematic features simulated in a realistic Arctic configuration between using the traditional and our suggested choices of EVP damping time scales, but it is minor considering model uncertainties associated with choices of many other parameters in sea-ice models

The Arctic Ocean in CMIP6 Models: Biases and Projected Changes in Temperature and Salinity

We examine the historical evolution and projected changes in the hydrography of the deep basin of the Arctic Ocean in 23 climate models participating in the Coupled Model Intercomparison Project Phase 6 (CMIP6). The comparison between historical simulations and observational climatology shows that the simulated Atlantic Water (AW) layer is too deep and thick in the majority of models, including the multi-model mean (MMM). Moreover, the halocline is too fresh in the MMM. Overall our findings indicate that there is no obvious improvement in the representation of the Arctic hydrography in CMIP6 compared to CMIP5. The climate change projections reveal that the sub-Arctic seas are outstanding warming hotspots, causing a strong warming trend in the Arctic AW layer. The MMM temperature increase averaged over the upper 700 m at the end of the 21st century is about 40% and 60% higher in the Arctic Ocean than the global mean in the SSP245 and SSP585 scenarios, respectively. Salinity in the upper few hundred meters is projected to decrease in the Arctic deep basin in the MMM. However, the spread in projected salinity changes is large and the tendency toward stronger halocline in the MMM is not simulated by all the models. The identified biases and projection uncertainties call for a concerted effort for major improvements of coupled climate models

Modified halocline water over the Laptev Sea continental margin : historical data analysis

Historical hydrographic data (1940s–2010) show a distinct cross-slope difference of the lower halocline water (LHW) over the Laptev Sea continental margins. Over the slope, the LHW is on average warmer and saltier by 0.2°C and 0.5 psu, respectively, relative to the off-slope LHW. The LHW temperature time series constructed from the on-slope historical records are related to the temperature of the Atlantic Water (AW) boundary current transporting warm water from the North Atlantic Ocean. In contrast, the on-slope LHW salinity is linked to the sea ice and wind forcing over the potential upstream source region in the Barents and northern Kara Seas, as also indicated by hydrodynamic model results. Over the Laptev Sea continental margin, saltier LHW favors weaker salinity stratification that, in turn, contributes to enhanced vertical mixing with underlying AW

Ocean Model Formulation Influences Transient Climate Response

The transient climate response (TCR) is 20% higher in the Alfred Wegener Institute Climate Model (AWI-CM) compared to the Max Planck Institute Earth System Model (MPI-ESM) whereas the equilibrium climate sensitivity (ECS) is by up to 10% higher in AWI-CM. These results are largely independent of the two considered model resolutions for each model. The two coupled CMIP6 models share the same atmosphere-land component ECHAM6.3 developed at the Max Planck Institute for Meteorology (MPI-M). However, ECHAM6.3 is coupled to two different ocean models, namely the MPIOM sea ice-ocean model developed at MPI-M and the FESOM sea ice-ocean model developed at the Alfred Wegener Institute, Helmholtz Centre for Polar and Marine Research (AWI). A reason for the different TCR is related to ocean heat uptake in response to greenhouse gas forcing. Specifically, AWI-CM simulations show stronger surface heating than MPI-ESM simulations while the latter accumulate more heat in the deeper ocean. The vertically integrated ocean heat content is increasing slower in AWI-CM model configurations compared to MPI-ESM model configurations in the high latitudes. Weaker vertical mixing in AWI-CM model configurations compared to MPI-ESM model configurations seems to be key for these differences. The strongest difference in vertical ocean mixing occurs inside the Weddell and Ross Gyres and the northern North Atlantic. Over the North Atlantic, these differences materialize in a lack of a warming hole in AWI-CM model configurations and the presence of a warming hole in MPI-ESM model configurations. All these differences occur largely independent of the considered model resolutions

Identifying climate change information needs for the himalayan region: Results from the GLACINDIA Stakeholder Workshop and Training Program

Here we present results of a workshop designed to bring together stakeholders from different states of the Indian side of the Himalayan arc and an international group of climate scientists in order to discuss how climate change research for this region can be tailored toward the needs of local communities. The stakeholder workshop was jointly organized by the Jawaharlal Nehru University (JNU), New Delhi, India, and the Climate Service Center 2.0, Hamburg, Germany, within the framework of the multidisciplinary international research project GLACINDIA. The project focuses on the water-related effects of changes in glacier mass balance and river runoff in western Himalayas. Given the research focus of the GLACINDIA project, the initial focus of the workshop was on glacier-related hydrological information. During stakeholder interactions the resulting discussion covered a much broader range of urgent climate change information needs for the Himalayan region.publishedVersio

Eddies and the distribution of Eddy Kinetic Energy in the Arctic Ocean

© The Author(s), 2022. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in von Appen, W.-J., Baumann, T. M., Janout, M., Koldunov, N., Lenn, Y.-D., Pickart, R. S., Scott, R. B., & Wang, Q. Eddies and the distribution of eddy kinetic energy in the Arctic Ocean. Oceanography, 35(2), (2022), https://doi.org/10.5670/oceanog.2022.122.Mesoscale eddies are important to many aspects of the dynamics of the Arctic Ocean. Among others, they maintain the halocline and interact with the Atlantic Water circumpolar boundary current through lateral eddy fluxes and shelf-basin exchanges. Mesoscale eddies are also important for transporting biological material and for modifying sea ice distribution. Here, we review what is known about eddies and their impacts in the Arctic Ocean in the context of rapid climate change. Eddy kinetic energy (EKE) is a proxy for mesoscale variability in the ocean due to eddies. We present the first quantification of EKE from moored observations across the entire Arctic Ocean and compare those results to output from an eddy resolving numerical model. We show that EKE is largest in the northern Nordic Seas/Fram Strait and it is also elevated along the shelf break of the Arctic Circumpolar Boundary Current, especially in the Beaufort Sea. In the central basins, EKE is 100–1,000 times lower. Generally, EKE is stronger when sea ice concentration is low versus times of dense ice cover. As sea ice declines, we anticipate that areas in the Arctic Ocean where conditions typical of the North Atlantic and North Pacific prevail will increase. We conclude that the future Arctic Ocean will feature more energetic mesoscale variability

On the role of horizontal resolution over the Tibetan Plateau in the REMO regional climate model

A number of studies have shown that added value is obtained by increasing the horizontal resolution of a regional climate model to capture additional fine-scale weather processes. However, the mechanisms leading to this added value are different over areas with complicated orographic features, such as the Tibetan Plateau (TP). To determine the role that horizontal resolution plays over the TP, a detailed comparison was made between the results from the REMO regional climate model at resolutions of 25 and 50 km for the period 1980–2007. The model was driven at the lateral boundaries by the European Centre for Medium-Range Weather Forecasts Interim Reanalysis data. The experiments differ only in representation of topography, all other land parameters (e.g., vegetation characteristics, soil texture) are the same. The results show that the high-resolution topography affects the regional air circulation near the ground surface around the edge of the TP, which leads to a redistribution of the transport of atmospheric water vapor, especially over the Brahmaputra and Irrawaddy valleys—the main water vapor paths for the southern TP—increasing the amount of atmospheric water vapor transported onto the TP by about 5. This, in turn, significantly decreases the temperature at 2 m by > 1.5 °C in winter in the high-resolution simulation of the southern TP. The impact of topography on the 2 m temperature over the TP is therefore by influencing the transport of atmospheric water vapor in the main water vapor paths. © 2018 Springer-Verlag GmbH Germany, part of Springer Natur